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DESIGN OF DRUGS INVOLVING 
RECEPTOR-LIGAND-DNA INTERACTIONS 

Cross-reference to Related Cases 

This application is a continuation-in-part of 
application serial number 08/158,689 filed on November 26, 1993 
which is hereby incorporated by reference in its entirety. 

Field of the Invention 

The present invention relates generally to rational 
drug design, in particular design of biologically active molecules 
using pharmacophores prepared according to the present 
invention. 

Background of the Invention 

Why certain chemical stractures and not others are 
present in nature has been a recurring question raised by scientists 
since the first organic natural products were characterized. Of 
equal interest has been elucidating what structural features within 
any given class of organic molecules are responsible for 
biological activity. Historically, the lack of satisfactory answers 
to both questions has relegated the development of biologically 
active molecules either to serendipity or to exhaustive synthesis 
and biological testing of large numbers of compounds. This 
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frastration is particularly evident in the pharmaceutical industry 
where the development of drug agonists and antagonists is often 
time consuming, tedious and expensive. 

This picture is beginning to change as more 
information is derived from modem molecular modeling 
techniques including characterization of the active sites in 
enzymes arid the ligand binding sites in receptors. Over the past 
15 years, another approach has emerged based upon a series of 
discoveries made with molecular models, wherein biologically 
active small molecules have been found to possess complementary 
stereochemical relationships with gene structure. This approach 
was first described in U.S. Patent No. 4,461,619 to Hendry, et al„ 
which is incorporated herein by reference. This simple molecular 
modeling technology was developed from observations first 
reported in 1977 of structural relationships between small 
molecules and nucleic acids, as described by Hendry, et al., / 
Steroid Biochem. Molec, Biol 42:659-670 (1992); Copland, et 
aL, /. Steroid Biochem. Molec. Biol 46:451-462 (1993); Hendry 
and Mahesh, /. Steroid Biochem. Molec. Biol 41:647-651 (1992); 
Witham and Hendry, J. Theor Biol 155:55-67 (1992); Hendry 
andMahesh, /. Steroid Biochem. Molec. Biol 39:133-146 (1991); 
Hendry, /. Steroid Biochem. 31:493-523 (1988); Lehner, et al., 
Molec. Endocrinol 1:377-387 (1987); Hendry, et al., /. Steroid 
Biochem. 24:843-852 (1986); Uberoi, et al.. Steroids 45:325-340 
(1985); Bransome, et al., /, Theor. Biol 112:97-108 (1985); 
Hendry, et aL, Proc. Natl Acad. ScL USA 78:7440-7444 (1981); 
and Hendry, et al., Perspect. Biol Med. 27:623-651 (1984), all of 
which are incorporated herein by reference. 
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The essential ingredient of all genes is a single, well 
defined polymer, deoxyribonucleic acid (DNA). DNA is a 
remarkably uncomplicated molecule composed of recurring 
sugar-phosphate units attached to one of four possible bases: 
adenine (A), thymine (T), cytosine (C) or guanine (G). The 
simplicity of gene structure is further evident in the Watson and 
Crick base pairing scheme of double-stranded DNA (A widi T 
and C with G), and the helical chirality (handedness) dictated by 
the absolute configuration of the sugar D-deoxyribose. Gene 
stmcture could conceivably be composed of many other chemical 
units, for example, other sugar stereoisomers such as L- 
deoxyribose or sugar homologs related to D-glucose. 

The products of gene structure, proteins, are also 
simple, ubiquitous molecules. Nature limits the structure of 
proteins by constmcting them from only twenty basic units, the 
amino acids; protein chirality is constrained by the absolute L- 
configuration of the amino acids. As in the case of nucleic acid 
subunits, a wide range of structural alternatives are possible for 
protein amino acids. Examples include changes in the chirality of 
a given aadno acid side chain (e.g., D-isoleucine), rearrangements 
in the pattern of atoms (e.g., the t-butyl isomer of isoleucine) or 
the addition of atoms (e.g., pipecolic acid, a homologue of 
proline). 

Structural constraints are also evident in the 
stereochemistry of low molecular weight natural products. 
Particularly conspicuous are limitations imposed by nature on the 
number, size, shape, elemental composition, and chirality of 
biologically active small molecules. For example, the pervasive 
neurotransmitters histamine and serotonin are unique in that 
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alternative structures with changes in the position or composition 
of heteroatoms and/or ring patterns generally do not exist in 
nature. Similarly, many small molecular weight hormones are 
few in number, have recurring structural patterns and possess a 
single absolute chirality. 

The source of the pervasive occurrence of 
physicochemical constraints on the structure of naturally 
occurring small molecules lies directly in the structure of die 
proteins which govern both their biosynthesis and bioactivity, i.e., 
enzymes and receptors, respectively. Ultimately, however, this 
stereochemical information is contained in the genes. According 
to the basic tenants of molecular biology, the information in DNA 
is replicated with remarkable precision and fidelity into newly 
synthesized DNA. It is also transcribed into RNA and 
subsequently translated into protein. 

This scenario, however, presents an apparent 
paradox. While the genetic template ultimately directs which 
proteins and small molecules are synthesized, as well as which 
proteins and small molecules will interact with each other, the 
undirectional flow of genetic information during translation 
suggests tiiat DNA stracture performs this function without 
recognizing the stracture of the small molecule. With few 
exceptions, such as certain antibiotics which bind directly to DNA 
and block transcription, small molecules are not considered to 
recognize or interact with the genetic template. Moreover, the 
structures of the molecules that are biosynthesized are thought to 
be unrelated to the stracture of the genes. 

In the initial search for stractural relationships 
between biologically active natural products and DNA, it became 
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apparent that the two-dimensional structures of DNA base pairs 
were analogous to many classes of small molecules, including 
gibberellic acid, a phytohormone; benzo [a] pyrene oxide, a 
carcinogen; the prostaglandin PGE2; morphine, a narcotic; 
estradiol, a hormone; riboflavin, vitamin Bi2; serotonin, a 
neurotransmitter; and actinomycin, an antibiotic. In addition to 
similarities in size and shape, numerous snuQl molecules contained 
donor/acceptor functional groups at locations where hydrogen 
bonds occurred between the base pairs. When overlaid on the 
base pairs, some compounds, such as the plant hormone 
gibberellic acid, the steroid hormone estradiol, and 
prostaglandins, contained heteroatoms separated by intemuclear 
distances similar to that of phosphate oxygens on adjacent strands 
of double-stranded DNA. This was particularly evident in 
functional groups attached at the 3 and 17p positions of the 
steroids. 

Using three dimensional Corey-Pauling-Koltun 
(CPK) space filling models, it became apparent that there were 
spaces between base pairs in partially unwound DNA that could 
accommodate a variety of small molecules. For example, 
estradiol could be inserted between base pairs in DNA, and the 
hydroxyl groups at 3 and 17p of estradiol were positioned such 
that they could form hydrogen bonds to phosphate oxygens on 
adjacent strands of DNA. Other steroids, including testosterone 
and progesterone, were also capable of stereochemical insertion 
between base pairs. In each case, complementary donor/acceptor 
linkages could be formed and the steroid conforaied well to the 
topography of the double helix. Attempts to insert any of the 
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non-naturally occurring steroid enantiomers into DNA resulted in 
poor fit in that donor/acceptor linkages were strained or could not 
form, and/or the overall sh^e of the molecules was incompatible 
with the helical topography of the DNA. 

5 Certain synthetic confounds with hormonal activity 

can also be accommodated within the DNA; in many cases, the fit 
of synthetic compounds such as diethylstilbestrol mimicked that of 
the natural hormone. In addition to mammalian steroids, 
prostaglandins, the insect hormone ecdysone and several 

10 phytohormones were also enable of stereochemical insertion and 

"recognition*' by the double helix. In the case of the plant 
hormone gibberellic acid, four stereospecific hydrogen bonds 
could be formed to donor/acceptor positions on the DNA, As 
with the steroids, only the naturally occurring enantiomer of 

15 gibberellic acid conformed to the topography of the double helix. 

One conclusion drawn from these studies is that 
certain chemical shapes, coupled widi heteroatom positioning 
compatible with that of the phosphate backbone of DNA and 
hydrogen bond positions of the base-pair template, potentiate 

20 partial or complete recognition between biologically active 

molecules and DNA. 

While it was possible to form complexes between 
DNA and a variety of molecules, amino acids did not initially 
show any clear accommodation to die space between base pairs. 

25 Certain compounds derived from amino acids, for example, 

neurotransmitters, fit into related sites. 

These relationships have been described as a 
stereochemical logic associated with gene structure. The 
stereochemical logic is defined as those unique features of nucleic 
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acid structure which ultimately dictate constraints on molecular 
structure, function, metabolism, and biologic activity. 

The use of molecular modeling as a tool to study 
organic structure has dramatically increased due to the advent of 

5 computer graphics. Not only is it possible to view molecules on 

computer screens in three dimensions but it is also feasible to 
examine the int^actions of ligands with various macromolecules 
such as enzymes and receptors, as reviewed by Borman, Chem. 
Eng. News 70:18-26 (1992). An abnost baffling array of 

10 software and hardware is now available and virmally all major 

pharmaceutical companies have computer modeling groups which 
are devoted to drug design. 

Modem methods of drug design include studies 
which focus on the binding of a molecule to a protein such as a 

15 polypeptide ligand for a receptor, or a steroid such as an estrogen 

or progesterone for a receptor. Similarly, drugs can be designed 
based upon the interaction of substrates with various enzymes. 
For the most part, however, binding sites in proteins have been 
difficult to characterize. There are many situations where other 

20 mechanisms must be involved to explain the feedback between 

protein regulation and regulation of gene expression. 

What is needed is a method for accurately predicting 
the biological activity of a given compound. The method should 
be easy to perform and should be able to predict both agonist and 

25 antagonist activity. 

Summary of the Invention 

The present invention is a method for identifying 
biological activity of molecules using pharmacophores. 
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According to the present invention, molecules are screened by 
determining the degree of "fit" in ±e pharmacophore. 

The method according to the present invention can be 
used to identify dmgs having increased biological activity or 
which have usefulness as antagonists or agonists, including, for 
example, estrogens and anti-estrogens. This method can also be 
used for the following: to predict the fit of compounds into 
nucleic acids, especially DNA; to predict the bioactivity of 
compounds, to screen compounds for toxicity; to design chemical 
groups to add to specific sites on molecules to facilitate 
metabolism or render the drag an agonist or antagonist; and to 
create molecules that mimic the activity of the DNA binding 
regions of receptors. 

The present invention also includes pharmacophores 
and the method of producing the pharmacophores and the use of 
the pharmacophores in predicting biological activity of a given 
compound. The present invention also includes the design of 
biologically active molecules using the pharmacophore. 

It is therefore an object of the present invention to 
provide a mediod which can be used to design a biologically 
active molecule. 

It is another object of the present invention to 
provide a method to screen and/or evaluate existing compounds 
for toxicological activity. 

Still another object of the present invention is to 
provide a method to predict the toxicity of compounds. 

Another object of the present invention is to provide 
a method to predict the toxicity of compounds for specific organs, 
tissues, and cells. 
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Yet another object of the present invention is to 
provide a method to design compounds that will have particular 
types of biological activities, including, but not limited to, 
hormonal, neurotransmitter, metabolic, genetic, immunologic, 
pathologic, toxic, and anti-mitotic activities. 

Still another object of the present invention is to 
provide a method to predict the bioactivity of compounds 
including, but not limited to estrogenic, anti-estrogenic, 
androgenic, anti-androgenic, progestational, anti-progestational, 
mineralocorticoid, retinoid, vitamin D like, thyroid, and 
glucocorticoid bioactivities. 

Yet another object of the present invention is to 
provide a method to create pharmacophores that can be used to 
design compounds such as drugs, hormones, neurotransmittCTS, 
agonists and antagonists more efficiently and economically. 

It is another object of the present invention to 
provide a method to create receptor pharmacophores that are 
molecular models of the portions of receptor molecules that bind 
to nucleic acids. 

Yet another object of the present invention is to 
provide a receptor phannacophore that can be used to design 
molecules that bind to nucleic acids with different affinity than the 
receptor. 

Yet another object of the present invention is to 
provide a pharmacophore that represents the three dimensional 
arrangement of solvent molecules around the ligand 
pharmacophore that binds to nucleic acids. 

Another object of the present invention is to provide 
a pharmacophore that represents the three dimensional 
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arrangement of solvent molecules around the receptor 
pharmacophore that binds to nucleic acids. 

It is yet another object of the present invention to 
provide a method to create a pharmacophore that is a three 
dimensional model of the nucleic acid binding domain of the 
receptor and of the ligand molecule that binds to the receptor and 
interacts with the nucleic acid at a different site. 

Another object of the present invention is to provide 
a phamiacophore that represents the three dimensional 
arrangement of molecules that can be attached to other 
pharmacophores to modify their biological activity. 

Still another object of the present invention is to 
provide a pharmacophore that represents the three dimensional 
arrangement of molecules that can be attached to other 
pharmacophores in order to design sites for enzymatic cleavage. 

These and other objects, features and advantages of 
the present invention will become apparent after a review of the 
following detailed description of the disclosed embodiments and 
the appended claims. 
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Brief Description of the Figures 

Figure 1 is a schematic of a cavity in DNA and the 
numbered sites which acconMnodate steroid/thyroid/vitamin A and 
D Ugands: testosterone (1,7); progesterone (2,7); aldosterone 
(2,5J and 9); Cortisol (2,3,4,5, and 7); estradiol (1 and 6); 
triiodothyronine (T3) (1 and 6); retinoic acid (6); and 1,25- (0H)2 
vitamin D3 (1 and 6); these are divided into two groups based on 
their interaction with either site 6 or 7. 

Figure 2 is a computer generated space filhng stereo 
view of the DNA cavity (A), which fits active estrogens oriented 
by energy calculations into the DNA cavity (B), whereas (C) 
shows the combined active surface of estrogens removed from the 
cavity in DNA that is used to construct the pharmacophore. 

Figure 3 demonstrates a volume contour map 
(yellow) in stereo with dummy atoms (magenta) surrounding the 
active molecules which were used in the construction of the 
pharmacophore (A); the empty pharmacophore (B); fit of the 
highly active estrogen 3,lip,17p-trihydroxy-7a-methylestra- 
l,3,5(10)-triene 11 -nitrate ester (hereinafter 7a-methylestradiol- 
liP-nitrate ester reported in Peters et al., 7. Med. Chem. 
32:2306-2310 (1989)) which is acconmiodated completely within 
the pharmacophore (C); and poor fit of the inactive estrogen 9p- 
estradiol which extends appreciably beyond the surface of the 
pharmacophore (D). 

Figure 4 shows a stereo view of volume maps (green 
mesh) and dummy atoms (magenta) for pharmacophores for 
estrogen (A), androgen (B), thyroid (C), anti-estrogen (D), and 
toxicity (E). Dummy atoms are not presented with die toxicity 
pharmacophore volume map. 
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Figure 5, left column, shows the volume maps (green 
mesh) and dummy atoms (magenta) for pharmacophores for 
estrogen (A), androgen (B), thyroid (C), anti-estrogen (D), and 
toxicity (E). The right colunm presents an orthogonal view (90 
degree rotation) of the corresponding images in the left column. 

Figure 6 demonstrates the relationship between the 
degree of fit of various compounds to the estrogen 
pharmacophore and the relative uterotropic (estrogenic) activity. 

Figure 7 shows the chemical structures of three 
separate molecules, SGI 100, SGI 101, and SGI 102 designed with 
pharmacophore technology. 

Appendix 1 is a QIC computer tape containing data 
files showing volume maps for each of the pharmacophores 
(estrogen, androgen, thyroid, antiestrogen, and toxicity) presented 
in Figures 4 and 5; the molecular database containing the dunmiy 
atoms for each pharmacophore designated by name; and a file 
containing the four color prints (Rgures 2-5) submitted with this 
application. 

Detailed Description of the Invention 

Many natural products have structures that exhibit 
stereochemical complementarity with nucleic acids, including 
amino acids, phytohormones, cycUc nucleotides, prostaglandins, 
insect hormones, steroid hormones, neurotransmitters, sugars, 
peptide hormones, thyroid hormones, pheromones, and vitamins. 
A striking example is the cyclopentanophenanthrene motif 
repeated in all classes of mammalian steroid hormones, such as 
estrogen and progesterone. Another example is the kaxurene 
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nucleus containing eight chiral centers which is evident in the 
gibberellin class of plant hormones. 

There are many ways to examine the stereochemistry 
of various configurations and conformations of nucleic acids. For 
example. Silastic polymer models can be constracted based upon 
computer derived space filling x-ray coordinates to reflect the 
stCTeochemistry of partially unwoimd DNA/RNA complexes, 
RNA-RNA complexes, ent-DNA (mirror image of DNA made 
with L-deoxyribose), and apurinic/apyrimidinic sites in DNA. 
These cavities reveal a sequence specificity in the fit of many 
molecules. The apurinic/apyrimidinic sites accommodate amino 
acids according to the known genetic code. The plant hormone 
gibberellic acid fits best into the partially imwound site 5'-dTdA- 
3', 5'-dTdA"3*; members of the mammalian steroid/thyroid 
hormone superfamily fit best into 5'dTdG-3', 5'-dCdA-3\ Each 
class of mammalian hormone forms xmique stereospecific donor- 
acceptor hydrogen bonds with DNA. The capacity to fit within 
these cavities in the manner of the index biologically active 
molecule correlates with the degree of biologic activity. It is not 
possible to fit chiral naturally occiuring molecules into ent-DNA. 
However, ent-DNA acconraiodates the biologic unknown chiral 
enantiomers, such as ent-progesterone. 

According to the present invention, computer 
modeling is used to examine the relationships between compounds 
and their fit in helical DNA. Although described herein with 
reference to double stranded, helical DNA, many of the same 
principles are applicable to double stranded RNA, and/or to RNA- 
DNA hybrids. Unless otherwise specified, double stranded RNA 
and DNA are to be considered equivalents as used herein. 
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Computer modeling can be used to view the interactions of 
molecules as well as to measure the energy of a given interaction. 
While a variety of software packages are available for computer 
modeling of molecules, a preferred software package is Sybyl 
software (version 6.03; Tripos Associates, St. Lx)uis, MO) for 
measuring the docking of various small molecular weight ligands 
into DNA. In the examples described herein, the software is run 
on a Silicon Graphics Indigo Extreme equipped with hardware 
stereo, i.e., Crystal Eyes (StereoGraphics, San Rafael, CA). 
Structures of small molecules are obtained via: the Cambridge 
Crystallographic Database, Lanfield Rd., Cambridge, England; 
construction with the Concord program or from fragment 
libraries and/or modifications of existing x-ray structures 
followed by energy nMuimization. All energy calculations are 
made using the Sybyl force field with a 1.2 A van der Waals 
parameter for hydrogen, as described by Hendry, et al., 7. Steroid 
Biochem. Molec, Biol 42:659-670 (1992) and Hendry, et al., /. 
Steroid Biochem. Molec, Biol 39:133-146 (1991). Charges are 
calculated using the Gasteiger-Huckel method to include a and % 
bonding. Partially unwound DNA cavities of various double 
stranded dinucleotide sequences are constructed from the Watson 
and Crick canonical B form of DNA by twisting each of the 
fourteen torsional angles along the sugar-phosphate backbone 
while maintaining the stractural integrity of hydrogen bonds 
between the base pairs. The Y-endo deoxyribose conformation of 
the sugars is employed and Kollman charges are calculated. 
Within thesexonstraints, attempts are made to adjust the degree of 
unwinding and the resulting degree of separation of the base pairs 
to accommodate the width of various ligands 
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Each ligand is insened into the cavity in DNA using 
van der Waals dot, mesh, and space filling surfaces in stereo to 
guide the docking procedure and minimize steric strain. The 
oxygens of the phosphate groups are permitted to act either as 
proton donors or acceptors and oriented to maximize the fit of 
any given ligand. The docking of the molecules is repeated 
several times. The distances between heteroatoms are monitored 
interactively to optimize the direction and distances of potential 
hydrogen bonds. While evaluating fit of compounds in a given 
hormone class, attempts are made to insert all candidate ligands 
into the DNA sequence with donor/acceptor linkages conmion to 
the hormone. Donor/acceptor relationships are further 
maximized by appropriate orientation of functional groups on the 
ligands, e.g., by adjusting the conformation of each structure to 
best n^mic the fit of the hormone. 

Van der Waals interactions of the candidate ligands 
are optimized with the DNA surface. The force field is used to 
calculate the relative fit of each ligand by assessing the optimum 
favorable energy change resulting from docking the ligand Each 
ligand is docked into the DNA and the steric fit is calculated from 
the change in van der Waals energy; the hydrogen bonding fit is 
calculated from the change in electrostatic energy using charges 
on donor hydrogens and acceptor heteroatoms. The energy 
changes are assessed for each ligand. The greater the negative 
energy change upon insertion of a given ligand into DNA, the 
more favorable the fit and the more stable the complex. Optimal 
docking is determined when no further increase in fit is observed. 
The most favorable change in energy is selected to compare the 
relative fit of each molecule. The overall fit of each ligand is 
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assessed by adding the change in kcal of the van der Waals and 
electrostatic energies and normalizing the fit to that of the parent 
hormone. 

Complementarity of Biologically Active Structures 

Although not wanting to bound by the following 
hypothesis, it is believed that the mechanism of action of 
hormones in the steroid/thyroid superfamily where ligand boimd 
receptor is known to regulate hormone responsive genes is 
mediated by physically inserting small molecules into DNA. The 
prior art teaches that the ligand causes a specij5c conformational 
change in the receptor protein which in turn contacts the DNA 
resulting in gene regulation. In contrast, according to the present 
invention, there is a receptor-mediated insertion of the ligand into 
DNA. This is consistent with the widely reported lack of 
correlation between hormonal activity and receptor binding for 
estrogenic steroids. According to Brooks, et al„ Recent Advances 
in Steroid Hormone Action, Moudgil, V.K. (ed) 443-466 
(Walther de Grayter, NY 1987), who made an extensive study of 
estrogen stracture-function relationships, receptor binding is 
essential for target-cell responses elicited by the steroids. 
However, the affinity of altered estrogens is not directiy related to 
the character or extent of the response. At the same time, the 
binding of estrogens and otiier steroids to DNA in the absence of 
receptor has been observed to be weak by numerous investigators. 
These include studies of the flat phytoestrogen coumestrol which 
might be expected to easily sUde between base pairs. In short, the 
bindmg of the steroid independently to either the receptor or the 
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DNA does not appear to be sufficient to fully explain hormonal 
response. 

Several pieces of evidence indicate that certain 
molecules, such as steroids, that have been shown to fit 
remarkably well between base pairs in DNA may elicit strong 
biological responses albeit through relatively weak interactions 
with DNA, Support for this concept is based on in vitro and in 
vivo experimental data and the energy calculations demonstrating 
a correlation between fit of estrogens into DNA and uterotropic 
activity, for example, studies of a new class of antitumor agents 
which have led to the discovery of drugs that are potent but act by 
binding weakly to DNA via intercalation, reported by Lee, et al., 
J. Med Chem. 35:258-266 (1992). Additional support that 
binding to receptors alone is insufficient to explain activity is 
provided by the observation of poor binding of a potent estrogen 
analog, ll^-acetoxyestradiol, which is considerably more active 
than estradiol. 

In general, given that degree of hormonal activity 
does not always correlate well with receptor binding but does 
correlate with fit into DNA, it follows that the mode of action of 
steroids involves a stage(s) in which the ligand is recognized by 
both the receptor and the DNA. This conclusion has also been 
reached in studies of anti-androgens. 

It is contemplated as part of the present invention that 
binding of the steroid to its receptor serves as a means to 
recognize the general class of hormone (e.g., estrogen versus 
androgen) whereas the type and degree of the fit of the steroid in 
the DNA is largely responsible for governing die magnitude of 
the biologic response. The steroid alone is incapable of proper 
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insertion into DNA either in vivo or in vitro without the receptor 
and thus would be unable to generate a full hormonal response 
without the receptor. A potentially critical role for the receptor 
upon binding to DNA, possibly in concert wi± other transcription 
factors, is to regulate the physicochemical properties of the site in 
DNA to permit insertion of the steroid, including the degree of 
unwinding, the capacity of heteroatoms to act as either donors or 
acceptors, and the pattern and degree of solvation. This is further 
supported by the decrease in the stnface hydrophobicity of the 
receptor upon binding estrogens and antiestrogens as well as 
phosphorylation which enhances binding of the estrogen receptor 
to specific DNA sequences, reported by Denton, et al., /. BioL 
Chem, 267:7263-7268 (1992). 

The process of a receptor-mediated insertion of 
steroids into DNA presumably involves several steps. For 
example, in the case of estrogens, the initial contact of the DNA 
by the steroid-receptor complex could involve a partially exposed 
D ring of the steroid with the A ring still attached to the receptor. 
Analysis of receptor binding data supports this possibility. 
Stereospecific recognition of the DNA by partial insertion and 
hydrogen bonding of the 17P-hydroxyl of the steroid D ring with 
the 5*-dTdG-3* strand could be followed by complete insertion 
and the recognition within the double helix manifest by the 
linkage of both stereospecific hydrogen bonds. The weak binding 
observed for the steroids with DNA suggests that the steroid/DNA 
complex might be short lived and rapidly reversible. Certain 
estrogen antagonists which are acconmiodated differently in DNA 
may form longer lived interactions within the site. 
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There are numerous possible scenarios and 
implications of receptor mediated binding of ligands to nucleic 
acids. For example, the lack of effectiveness of certain anti- 
estrogens in tissues which lack estrogen receptor might be due in 
part by die inability of the antagonist to be transported to the 
DNA obviating the insertion of the ligand. Mutation of the site 
which accommodates the estrogens would also result in improper 
recognition of the ligand and would be predicted to no longer 
permit the gene to be properly regulated either by agonists or 
antagonists. There might be multiple levels at which a given 
compound might act, as suggested by the observation that cavities 
in double stranded RNA and RNA-DNA hybrids can 
accommodate various ligands, for example, estradiol in 5'-rUrG- 
3*,5'-dCdA-3\ 

Molecular Modeling 

Molecular modeling was performed with 
Sybyl/Mendyl 5.4 (Tripos Associates, St. Louis, MO) using an 
Evans and Sutherland PS390 graphics computer equipped with a 
stereographic viewer. Stmctures of piperidinedione ligands were 
provided via constmction with the Concord program or from 
fragment libraries followed by energy minimization. Energy 
calculations were made with Sybyl/Mendyl force field and a 1 .2A 
van der Waals parameter for hydrogen. Charges were calculated 
using the Gasteigner-Huckel method which includes a and tc 
bonding. Partially unwotmd DNA was constmcted witii 3*-endo 
deoxyribose conformation and Kollman charges. This method has 
been described by Hendry, et al., 7. Steroid Biochem. Molec. Biol 
42:659-670 (1992); Hendry, et aL, 7. Steroid Biochem. Molec. 
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BioL 39:133-146 (1991), and Hendry et al., / Steroid, Biochem, 
Molec. BioL 49: No. 4-6, pp. 269-280 (1994) the teachings of 
which are hereby incorporated by reference in their entirety. 

The ligands were inserted into the cavity in DNA 
using van der Waals dot surfaces and the stereoviewer to guide 
the docking procedure and minumze any obvious steric strain. 
The distances between heteroatoms were monitored interactively 
to optimize the direction and distances of potential hydrogen 
bonds. Donor/acceptor relationships were further maximized by 
appropriate orientation of functional groups on the ligands, e.g., 
by adjusting the conformation of each stmcture. Attempts were 
made to optimize van der Waals interactions of the candidate 
ligands widi the DNA surfaces. The force field was used to assess 
the relative fit of each ligand by quantitating the optimum 
favorable energy change resulting from docking the ligand. 
Steric fit was calculated from the change in van der Waals energy, 
the hydrogen bonding fit was calculated from the change in 
electrostatic energy using charges on donor hydrogens and 
acceptor oxygens. The greater the negative energy change upon 
insertion of a given ligand into DNA, the more favorable the fit 
and the more stable the complex. Docking was completed when 
no further increase in fit was observed. The most favorable 
change in energy was selected to compare the relative fit of each 
molecule. The overall fit of each ligand was assessed by adding 
the change in kcal of the van der Waals and electrostatic energies 
and normalizing the value to that of the best fitting molecule 
(100%). It should be noted that while the energies reported here 
were derived from widely used force field calculations, they were 
not en5)irically derived. Thus, the absolute values in kcal do not 
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have independent experimental significance. At the same time, 
they are valuable indicators of the relative degree of fit into DNA 
of candidate molecules. 

Previous studies using space filling models indicated 
that 3-phenylacetylamino-2,6 piperidinedione was capable of fully 
inserting between base pairs in DNA and forming a stereospecific 
hydrogen bond between the imino proton of the piperidinedione 
ring and a negatively charged phosphate oxygen on the 
deoxyribose-phosphate backbone. Results employing computer 
graphics confirmed this observation. Energy calculations further 
demonstrated that this compound had favorable van der Waals 
contacts of approximately -17.7 kcal when inserted into DNA 
with an electrostatic energy of approximately -21.7 kccd resulting 
from the stereospecific hydrogen bond (2.7 A) to phosphate. 
Increased fit of the ligand was obtained by substituting a para 
hydroxl group on the phenyl ring; this substitution enabled a 
second hydrogen bond to be formed between the hydroxyl group 
and a phosphate oxygen on the adjacent DNA strand. The 
increase in fit measured by energy calculations due to the second 
hydrogen bond (2.64 A) was reflected in an additional -24.6 kcal 
in electrostatic energy. Other substitutions which were made on 
the 3-phenylacetylamino-2, 6-piperidinedione skeleton did not 
significantly increase fit demonstrated by the normalized energy 
calculations for certain halogenated analogs. 

Synthesis 

The synthesis of the unsubstituted derivative, 3-N- 
phenylacetylamino-2, 6-piperidinedione has been briefly 
described by Burzynski, et al.. Drugs of the Future 10:103 



SOBSTITUTE SHEET (RULE 26) 



PCTAJS94/13765 

22 

(1985), and was used as a general method for the preparation of 
the desired compounds. Appropriate phenylacetic acids were 
reacted with N-hydroxysuccinimide in the presence of NJ^- 
dicyclohexylcarbodiimide (DCC) which gave succinimide esters. 
The active esters were stable enough to isolate for physical and 
spectroscopic characterization although the major portions of the 
esters were used for the next reaction without purification. The 
active esters were reacted with L-glutamine in the presence of 
sodium bicarbonate to obtain the glutamine derivatives. However, 
due to the difficulties of obtaining the analytical samples, crude 
products were directly used for the next reaction. To prepare 
active esters, the glutamine derivatives were again reacted with N- 
hydroxysuccinimide in the presence of DCC to give the active 
esters which without purification were heated at 95-100"C to 
obtain die desired 2,6-piperidinediones in various yields. During 
the heating process the compoimds were racemized. 

Biological Evaluation 

These synthetic derivatives were assessed for 
biological potency by measuring their growth inhibitory effects 
on various cell lines using concentrations of 4 nM based on the 
reported IC50 of 3-phenylacetylamino-2,6 piperidinedione in Nb2 
cells. In YAK lymphoma cells, the p-hydroxy compound was the 
most active derivative. This compound, p-hydroxy-3- 
phenylacetylamino-2,6-piperidinedione was also the most active 
analog when tested in human leukemia (K652) cells. A dose 
response comparison in K562 cells showed that it was more active 
than the unsubstituted compound over the concentration range 
tested (10"5 to 10'^ M). Prolactin stimulated growth of rat Nb2 
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lymphoma cells was inhibited by each of the compounds with p- 
hydroxy-3-phenylacetylamino-2,6-piperidinedione manifesting the 
greatest activity. Compound p-hydroxy-3-phenylacetylamino-2,6- 
piperidinedione was more active in Nb2 lymphoma cells than 3- 
phenylacetylamino-2,6 piperidinedione over the range tested (lO'^ 
M to 10-3 M). 

Further analysis of growth inhibition of the most 
active analog p-hydroxy-3-phenylacetylamino-2,6-piperidinedione 
compared with the parent compound p-hydroxy-3- 
phenylacetylamino-2,6-piperidinedione was performed in MCF-7 
(E-3) human breast cancer cells. Both 3-phenylacetylamino-2,6 
piperidinedione and p-hydroxy-3-phenylacetylaniino-2,6- 
piperidinedione inhibited estrogen simulated cell growth. In a 9 
day model, p-hydroxy-3-phenylacetylamino-2,6-piperidinedione 
was more active than 3-phenylacetylamino-2,6 piperidinedione 
with IC50 comparable to tamoxifen (3-phenylacetylaimno-2,6 
piperidinedione, 3x10'^ M; p-hydroxy-3-phenylacetylamino-2,6- 
piperidinedione, 7xl0'6 M, tamoxifen, lxl0"7 M). The open 
chain hydrolysis product of 3-phenylacetylamino-2,6 
piperidinedione, PAG, did not inhibit cell growth even at high 

concentrations (i.e., lO'^M). 

Computer modeling coupled with energy calculations 
confirm that 3-phenylacetylamino-2,6-piperidinedione is capable 
of inserting between base pairs in partially unwound double 
stranded DNA and forming an energetically favorable complex. 
A hydroxyl group placed in the para position of the phenyl ring 
of 3-phenylacetylamino-2,6-piperidinedione enabled formation of 
a second hydrogen bond thereby linking both DNA strands. This 
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added hydrogen bond resulted in a greater fit in the DNA as 
assessed by energy calculations, i.e., 3-phenylacetylannno-2,6- 
piperidinedione (61%) versus the p-hydroxy derivative p- 
hydroxy-3-phenylacetylainino-2,6"piperidinedione (1 00%). 
Various substitutions at the para position as well as ortho and 
meta positions with fluorine and chlorine did not result in a 
significant increase in fit compared to 3-phenylacetylamino-2,6- 
piperidinedione. 

When the analogs of 3-phenylacetylamino-2,6- 
piperidinedione were synthesized and examined for the capacity to 
inhibit cancer cell growth, the p-hydroxy derivative, p-hydroxy- 
3-phenylacetylainino-2,6-piperidinedione, which was predicted to 
be the most active compound based upon fit into DNA was 
consistently found to be the most potent compound. That die 
capacity of the hydroxy group of p-hydroxy-3- 
phenylacetylamino-2,6-piperidinedione to form a second 
hydrogen bond to DNA was responsible for the predicted increase 
in activity is further supported by the lack of increased potency of 
p-fluoro derivative as well as other halogenated derivatives which 
were incapable of forming analogous hydrogen bonds. These 
observations demonstrate that among the compounds examined, a 
conelation exists between degree of fit into DNA and predicted 
biological potency. These findings also support the contention 
that stereochemical complementarity of small molecules widi 
nucleic acids can be a powerful tool for designing new drugs. 

The mode of action of these piperidinediones is still 
not proven. One mode of action might involve insertion into 
DNA as suggested by the computer modeUng results and the 
observation that DNA synthesis measured by thymidine 
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incorporation was significantly inhibited upon treatment with 3- 
phenylacetylamino-2,6-piperidinedione in Nb2 lymphoma cells. 
No covalent adducts of 3-phenylacetylamino-2,6-piperidinedione 
with DNA were detected and the binding was observed to be weak 
and reversible, in con^arison to classical intercalating drugs. 
Another possible mode of action is suggested by the observation 
that compound p-hydroxy-3-phenylacetylamino-2,6- 
piperidinedione inhibits estrogen stimulated cell growth in MCF-7 
ceils, which is comparable to that of the established anti-estrogen 
tamoxifen. In contrast to tamoxifen, however, neither 3- 
phenylacetylamino-2,6'piperidinedione nor p-hydroxy-3- 
phenylacetylamino-2,6-piperidinedione exhibited appreciable 
binding for the estrogen receptor. At the same time, direct 
binding of such anti-estrogens to DNA appears to be weak. Taken 
as a whole, these findings support the possibility of a weak 
interaction of the piperidinediones with both DNA and the 
estrogen receptor, involving a receptor mediated insertion of the 
ligand into DNA. 

Development of Pharmacophores 

Molecular modeUng, as described above, faciUtates 
the establishment of the best fit of molecules into nucleic adds 
such as double-stranded DNA based on steric and electrostatic 
considerations. Individual molecules, such as estradiol, fit 
optimally into specific sites on DNA based on the location of 
specific nucleotides and the bonding characteristics of individual 
heteroatoms (see Example 2). Molecules that are related to a 
specific molecule such as estradiol but display chemical 
differences will fit into the estradiol site with different degrees of 
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precision: some may fit better and give rise to estradiol agonistic 
responses while those with poor fit display weak estrogenic 
activity. These different molecules may be aligned relative to the 
docking of heteroatoms with heteroatoms on the DNA to optimize 
electrostatic interactions. In the creation of pharnfiacophores 
described below, molecules with activity equal to or greater than 
that of the hormone are chosen for alignment. To date, such 
molecules fit equally well or better than the hormone into DNA 
using the energy calculation methodology described above. 
Molecules which do not fit as well into DNA as the parent 
hormone are excluded from inclusion in the constraction of the 
pharmacophore. The alignment of the combined surfaces of the 
molecules occupies a specific volume of space thereby forming a 
three dimensional shape. 

Phannacophores are three dinaensional arrangements 
of chemical groups related to a given biological activity which 
enables meaningful comparison of molecules exhibiting the same 
biological function (Naruto et al., Eur. /. Med. Chem. 20:529-532 
(1985)). Pharmacc^hores can be derived by simple overlap of 
active stractures or common functional groups in the molecules. 
Without a way to orient the molecules e.g., based upon fit with 
another macromolecule - a receptor, enzyme, or in this case 
DNA, it is difficult and, in some cases, impossible to construct a 
reliable pharmacophore. This problem results in part from the 
fact that even closely related active molecules frequently fit into 
macromolecules in very different ways. 

A pharmacophore, as used herein, is defined as a 3- 
dimensional shape having a specific volume derived from the 
combined van der Waals surface of active molecules oriented by 
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fit into DNA, coupled with point charges located adjacent to the 
surface. A pharmacophore represents an aggregate array of 
positions in space of a series of molecules having the same or 
similar biological activity. The van der Waals surface can be 
represented in various ways including as a volume map, a dot 
surface, or a Connolly surface. The point charges are represented 
as dummy atoms whose positions are determined by the average 
positions of functional groups on active molecules which can form 
hydrogen bonds. Suitable charges are placed on the dmmny atoms 
consistent with the capacity of the active molecules to form 
hydrogen bonds. The pharmacophores are specific for different 
conq)oimds, their related molecules and a particular biological 
activity. According to the present invention, within the general 
class of molecules called hormones, an estrogen pharaaacophore, 
an anti-estrogen pharmacophore, an androgen pharmacophore, a 
thyroid hormone pharmacophore, and a toxicity pharmacophore 
(shown in Figures 2-5 and Appendix 1) have been disclosed. It 
should be emphasized that these created pharmacophores do not 
exist as such in nature and are the product of aligning several 
related molecules to conmion binding sites in DNA using methods 
as described herein. 

Many other pharmacophores have been constracted 
using the method described in this application. A pharmacophore, 
once created, stands alone and is subsequentiy independent from 
the nucleic acid that was involved in its formation. Thus, after 
formation of the pharmacophore, one no longer needs to use the 
DNA as a template for the design of biologically active molecules. 
The pharmacophore itself can be used to generate new molecules 
that will possess the same or similar stractural and charge features 
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that are represented by the pharmacophore. This is a completely 
different concept from the one of using the DNA as the model for 
the design of compoimds. The pharmacophore can be used itself 
for any number of apphcations, including but not limited to the 
following: as a screening tool for drug development; to determine 
if a particular compound will possess bioactivity of a certain type, 
for instance estrogenic or androgenic activity; for toxicological 
evaluation; and to design compounds that possess increased or 
decreased binding affinity for DNA. 

Each pharmacophore has a characteristic shape, 
topology, volume, and electrostatic profile, A pharmacophore is 
accurately described by its three dimensional shape which is 
represented by a coordinate system that is configured in computer 
memory (see Figures 2-5 for examples of pharmacophores which 
are also contained in data files in Appendix 1, the magnetic QIC 
tape submitted with the application). Each specific atom within a 
molecule that fits in a pharmacophore has a specific location 
relative to the docking heteroatoms. The individual atoms also 
have electrical charges assigned to them. These charges are 
represented nimierically and through many other ways including 
the use of colors and shading to indicate field strength. As the 
degree of steric and electrostatic fit between the pharmacophore 
and the dummy atoms increases, resulting in a negative energy of 
interaction (-kcal), the efficacy of the pharmacophore increases 
which could manifest as increased bioactivity. The term "energy 
of interaction" as used herein is the total energy in -kcal of a 
molecule as it is being fitted into a pharmacophore. This has been 
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observed in the case of molecules that fit within the estrogen 
pharmacophore and bioactivity in a uterotropic assay. The 
volume of a pharmacophore is described in cubic angstroms. The 
pharmacophore can be cross sectioned precisely in any plane and 
internal distances measured with an angstrom ruler. The 
circumference of any cross section is easily measured with 
morphometric analysis. Similarly, specific subregions of the 
pharmacophore, such as tihe site that binds to the DNA, can be 
subjected to the same methods of analysis. 

Construction and Utility of Pharmacophores 

An example of ±e construction and utilization of the 
estrogen pharmacophore is given below. Hgure 2A is a computer 
generated space filling stereo view of the DNA cavity which fits 
estrogens. The fit of active estrogens oriented by energy 
calculations into the DNA cavity, using the methods described 
above, is presented in Figure 2B. Figure 2C shows the combined 
active surface of estrogens removed firom the cavity in DNA that 
is used to construct the pharmacophore. The atoms are colored in 
the following manner: carbon/white; hydrogen/cyan; 
nitrogen/blue; oxygen^ed; phosphoras/yellow. Figure 3 
demonstrates a volume contour map (yellow) in stereo with 
dummy atoms (magenta) surrounding the active molecules which 
were used in the construction of the phannacophore (A); the 
empty pharmacophore (B); fit of the highly active estrogen 
3,1 ip,17p-Trihydroxy-7a-methylestra-l,3,5(10)-triene 1 l-nitrate 
ester (hereinafter 7a-methylestradiol-lip-nitrate ester reported 
in Peters et al., 7. Med Chem. 32:2306-2310 (1989)) which is 
acconmiodated completely within the pharmacophore (C); poor fit 
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of the inactive estrogen 9p-estradiol which extends appreciably 
beyond the surface of the pharmacophore (D). Rgures 4 and 5 
present examples of the three dimensional appearance of estrogen, 
androgen, thyroid, anti-estrogen, and toxicity pharmacophores. 
The data files for the volmne maps and dmnmy atoms for each of 
these pharmacophores and for Figinjes 2 - 5 are contained in 
Appendix 1, the magnetic QIC 120 tape submitted with the 
application. This QIC computer tape has the following format and 
files generated from Sybyl 6.0 run on a Silicon Graphics Indigo 
Extreme. 

1 ) /usr/people/guest/empty .dir/estrogCTnpharmacophore2.dsp 
This is the volume map for the estrogen pharmacophore. 

2) /usr/people/guest/empty.dir/androgenpharmacophore,dsp 
This is the volume map for the androgen pharmacophore, 

3) /usr/people/guest/empty.dir/thyroid pharmacophore.dsp 
This is the volume map for the thyroid pharmacophore. 

4) /usr/people/guest/empty.dir/antiestrogenpharmacophore.dsp 
This is the volume map for the antiestrogen pharmacophore. 

5) /usr/people/guest/empty.dir/toxicitypharmacophore.dsp 
This is the volume map for the toxicity pharmacophore. 

6) /usr/people/guest/en5)ty.dir/dummyatoms.mdb 

This is the molecular database file containing the appropriate 
dummy atoms for each pharmacophore designated by name. 

7) /usr/people/guest/empty.dir/patent94final 

This is a file containing the four color prints (Figiu-es 2 through 
5) submitted with the patent application describing die 
pharmacophores constructed with the Silicon Graphics program 
Showcase. Copies of the Appendix 1 tape containing the data files 
are also retained by the inventor. These files contain die 
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information ±at represents the three dimensional configuration of 
each pharmacophore. These files may preferably be read with 
Sybyl software (version 6.03; Tripos Associates, St. Louis, MO) 
running on a computer system such as; CONVEX, ESV, IBM 
5 RISa6000, IRIS 4D, SUN 4, or VAXA^S, equipped with a QIC 

t^ drive. In addition, other media such as 4 and 8 mm DAT 
tapes may be used to store and read these data files on the systems 
mentioned above. 

Quantitative measurements of the degree of fit of 

0 various con:q)ounds to the pharmacophore are shown in Hgure 6. 
Fit is determined by measuring the amount of volume of each 
structure which could be placed widiin the pharmacophore 
volume map and normalizing the value to that of the natural 
hormone estradiol set at 50%. Electrostatic interactions with the 

5 dummy atoms are optimized for each compound and calculated 

using tiie Tripos force field. The electrostatic energy value is 
normalized to that of estradiol set at 50%. In this study, the total 
of volume fit and electrostatic fit are treated equally and totaled to 
reflect the overall fit in the pharmacophore. As shown in Figure 

) 6, degree of fit to the pharmacophore correlates highly to relative 

uterotropic (estrogenic) activity. In comparison to estradiol, 7a- 
mediylestradiol-llp-nittate ester (labeled 2), which is not part of 
the data set used to constmct die pharmacophore, fits appreciably 
better than estradiol Qabeled 1). In contrast, 7a-methylestradiol- 

1 1 ip-nitrate ester binds poorly to the estrogen receptor (less than 
6% of the binding of estradiol). The uterotropic values for 7a- 
metiiylestiradiol-lip-nitrate ester (labeled 2) relative to estradiol 
set to a normalized value of 100% are considerably greater than 
tiiat of estradiol (Peters et al., /. Med. Chem. 32:2306-2310 
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(1989)) as predicted by fit into the phannacophore. In contrast, it 
is not possible to fit 9p-estradiol (labeled 3) into the 
phannacophore and as predicted this analog has little uterotropic 
activity. In summary, fit of compounds to the pharmacophore 

5 correlates with biological activity and can thus be used to design 

new compoimds by virtue of their fit. It is noteworthy that the 
highly potent estrogenic nitrate ester binds very poorly to the 
estrogen receptor Thus, it would not be possible to predict the 
estrogenic activity of this analog on the basis of receptor binding 

10 or from a pharmacophore derived from a putative binding site in 

a protein receptor. In many cases, compounds with greater 
estrogenic activity than the natural hormone estradiol bind 
relatively poorly to the estrogen receptor. 

The fit of compounds into DNA is consistent with, 

15 but not the same as, fit into the pharmacophores. The distinction 

between the DNA cavities and the pharmacophores is that the 
surfaces were derived from different structures, i.e. the DNA 
cavities from the DNA structure and the pharmacophores from 
the combined surfaces of active compounds. As shown in the 

20 examples and in Figures 2-5, and in datafiles contained on the 

magnetic tape as Appendix 1, the degree of fit to the 
pharmacophores is obtained by fit to the three dimensional map 
which represents the active compounds. The pharmacophores 
enable quantitative determination of the degree of fit to the 

25 combined surfaces of the active compounds and this information 

cannot be obtained from fit into DNA. Moreover, the fit of 
active compounds to the pharmacophore can be quantitated based 
upon the portion of the molecule which does not fit within the 
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pharmacophore volume. This enables automatic assessment of 
predicted inactive structures. 

Toxicity Pharmacophores (Toxicophores) 

Pharmacophores can be constructed to represent a 
three dimensional sh^e that is predictive of toxic biological 
activity. Such pharmacophores, called toxicophores, have regions 
that would potentially damage DNA. A typical toxicophore has 
been constructed using, tetrodotoxin, dioxin, RU 486, dilantin, 
thaUdomide and oroflex, among other compotmds. An example of 
this toxicophore is provided in Figures 4E and 5E and on the 
submitted magnetic data tape as Appendix 1. For example, by 
overlaying this toxicophore on another pharmacophore such as 
the estrogen pharmacophore (Figure 4A and 5A), a drug designer 
would know to avoid designing an estrogenic compound with 
certain molecular groups that might impart toxic activity if these 
groups extended into the three dimensional space occupied by the 
toxicophore. This approach would greatly facilitate and 
economize drug design by guiding the designer to avoid 
synthesizing estrogenic compounds that might have damaging 
effects on DNA as opposed to proceeding with syndiesis and 
purification and subsequently discovering that the compound 
possesses dangerous toxicity. 

Solvent Pharmacophores (Aquaphores) 

Pharmacophores in their relationship to nucleic acids 
are usually surrounded by a solvent. The predominant solvent in 
living organisms is water and accordingly, most pharmacophores 
exist in an aqueous environment. Water is the preferred 
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embodiment of the solvent pharmacophore and is termed an 
aquaphore. Pharmacophores, and their molecules may also be 
placed in non-aqueous environments for various purposes such as 
crystallographic studies or other analytical procedures. 
5 In Uving organisms, the aqueous environment 

surrounding the pharmacophore also has an intimate association 
with the adjacent nucleic acid. This aqueous shell assists in the 
optimal fit of the pharmacophore into the cavity of the double 
stranded DNA, and has its own three dimensional shape. The 

10 optimal steric and electrostatic placement of water molecules in 

the space between the pharmacophore and the DNA is achieved in 
the present invention. This three dimensional shape is called a 
solvent pharmacophore, and can be described in all the ways listed 
above for the pharmacophores based on other molecules such as 

15 estrogen. Solvent pharmacophores assist the designer of 

confounds by placing limits on the dimensions of a compoxmd 
designed using a particular pharmacophore as a template. In 
addition, solvent pharmacophores assist the creator of 
pharmacophores because the solvent shell or cage represented by 

20 the pharmacophore provides enhanced ability to properly align 

molecules relative to DNA during the creation of the 
pharmacophore. 

Receptophore 

25 Many molecules, such as steroid hormones, are 

shuttled to the nucleus by other molecules known as receptors 
asai and O^Malley, Ann. Rev. Biochem, 63:451-486 (1994)). 
These receptors bind the hormones (called ligands), bind to the 
nucleic acids, for example in their DNA binding domain, and 
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present ligands to nucleic acids such as DNA. EAadence suggests 
that the binding of the receptor to the DNA causes a 
conformational change in the DNA to facilitate insertion of the 
ligand (NardulU et al., Molec. Endocr. 7:331-340 (1993)). The 
pharmacophore concept is based on the three dimensional shape of 
the optimal fit of related molecules into nucleic acids such as 
partially unwound, double-stranded DNA. The DNA binding 
domain of the receptor can be modeled into a three dimensional 
shape based on the same principles described above for the 
pharmacophore. The resultant shape is termed a receptophore and 
is the three dimensional representation of the sites of interaction 
of the receptor and the nucleotides of the DNA. The DNA binding 
region of each receptor likely gives rise to a different 
receptophore. This receptophore provides a valuable tool to 
molecular designers interested in developing new receptors, or in 
modulating receptor binding to DNA. 

Receptophore-Pharmacophore Pairs 

The nucleic acid binding region of receptors and 
their ligands can be modeled as * receptophores and 
phamiacophores, respectively. The configiuration of the 
receptophore and its associated pharmacophore in their proper 
ahgnment relative to their respective DNA binding regions 
constitutes a specific pair of shapes that represents the minimal 
molecular unit for DNA binding and ligand insertion. Designers 
of compounds utilize this information to synthesize and screen 
molecules to modify the facility of docking and ligand insertion. 
Such modifications may provide a host of new therapies such as 
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treatments for hormone dependent carcinomas of the prostate or 
breast. 

Metabophores 

Most naturally occurring confounds are derived 
from antecedents or precursors in a synthetic pathway and also 
are destined for inactivation in a catabolic pathway. Many 
precursors and metabolites of confounds are less active due to 
the addition of an extra group such as a methyl group or 
acetylation of a specific site. In some cases, precursors and 
metabolites of a molecule have different groups added 
sequentially to the same site on the active molecule, creating a side 
chain. Knowledge of a site of preferred addition or deletion of 
chemical groups assists in the design of molecule with enhanced 
or reduced activity. 

These sites can be modeled relative to the 
phannacophore to produce a three dimensional representation of 
a preferred site for modification of the molecule. This three 
dimensional representation, termed a metabophore, provides 
constraints for rational design of active and inactive variants of 
the parent molecules that fit into the pharmacophore. Analysis of 
which chemicals can effectively be added at the attachment point 
of the metabophore to the pharmacophore reveals the most 
favorable molecules to pursue for syndiesis, purification and 
testing. 

It will be appreciated that other embodiments and 
uses will be apparent to those skilled in the art and that the 
invention is not limited to these specific illustrative examples. 
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Example 1 

Fit into partially unwound double stranded DNA 
using the mammalian female hormone estradiol 

Computer modeling has demonstrated that the 
mammalian steroid progesterone is a remarkable "lock and key" 
fit into DNA at 5'-dTdG-3\ 5^-dCdA-3' (Figure 1). Each of the 
known x-ray crystal stractures of progesterone is capable of 
forming two stereospecific hydrogen bonds and a stable complex 
measured by force-field calculations. Remarkable 
complementarity is evident in the complex by the overlap of 
hydrophilic and hydrophobic regions of the steroid and DNA* 
The enantiomer of progesterone which does not occur in nature 
does not fit. The plant hormone gibberellic acid has also been 
shown to fit between base paire but in a different sequence, i.e., 
5'-TdA-3\ 5'-dTdA-3\ Four stereospecific hydrogen bonds are 
formed within the couples: ent-gibbCTellic acid does not fit. 

The mammalian female hormone estradiol also fits in 
DNA (Figure 1). Two stereospecific hydrogen bonds of 
approximately 2.65 A are formed between each hydroxyl group 
of the steroid and phosphate oxygens on adjacent strands. The 
overall fit within the complex is about -59 kcal. Most alterations 
of the positions of the hydroxyl groups on the estratriene nucleus 
result in a substantial loss of potential electrostatic interactions 
with the DNA. Moreover, most alterations of the absolute 
stereochemistry of the cyclopentanophenanthrene ring pattem also 
result in a substantial loss of potential electrostatic interactions 
with the DNA. Moreover, most alterations of the absolute 
stereochemistry of the cyclopentanophenanthrene ring pattem also 
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result in a poor fitting molecule. This is evident in attempting to 
fit 9P-estradiol into DNA. The puckering in the steroid caused by 
inversion of the estradiol stereochemistry at C-9 from a to p 
prevents complete insertion between the base pairs. Even if 
forced into DNA witiiout regard to strain caused by the overlap 
of van der Waals surfaces, 9P-estradiol can form only a single 
hydrogen bond. The relative fit into DNA resulting from partial 
insertion of 9p-estradiol is about -17 kcal. 9p-estradiol is inactive 
when tested in vivo for estrogenic (uterotropic) activity. 

The finding that estradiol is a "lock and key" fit into 
DNA, although most stmctural alternatives to estradiol fit poorly 
provides further support for the premise that DNA 
stereochemistry contains the master blueprint for natural product 
structures. 

Example 2 

Correlation of Estrogenic Activity with the Fit of 
Estrogens and Related Analogs into DNA. 

That fit into DNA measured by energy calculations 
can be correlated with biologic activity was demonstrated using a 
series of estrogens and related sjmflietic analogs. The molecules 
which are inactive in uterotropic assays fit poorly into DNA. 
Molecules diat fit into DNA better than estradiol are all more 
active than estradiol in the uterotropic assays, for example, llp- 
acetoxyestradiol (approximately - 68 kcal). This correlation is 
also observed with nonsteroidal, synthetic estrogens, such as the 
potent synthetic estrogen, transdiediylstilbestrol (approximately - 
62 kcal), which fits well, whereas the poorly active geometrical 
isomer cis-diethylstilbestrol (approximately -20 kcal) is a poor fit. 
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Example 3 

Biosynthetic Pathways Reflect Increasing Fit into 
DMA Whereas Jnactivation Pathways Lead to Decreasing Fit into 
5 DNA. 

Molecular modeling studies conducted with the 
mammalian hormone progesterone and the plant hormone 
gibberellic acid have shown that each step in the respective 
bios5rathetic pathway reflects a stmctural change that results in 

10 increased fit in DNA. For example, when considering possible 

stereoisomers which could result in any given step in 
progesterone biosynthesis, the best fitting structure is one which 
was produced in nature. In sharp contrast, each step in the 
inactivation of progesterone eventually leading to the excreted 

15 glucuronides and sulfates resulted in the worst possible fitting 

stereoisomers. 

The two possible dihdro reduced metaboUtes of the 
male hormone testosterone were examined for fit and correlation 
with biological activity. Relative to testosterone (100%), 5p- 

20 dihydrotestosterone is a poor fit (84%), whereas its epimer 5a- 

dihydrotestosterone fits even better than testosterone (102%). 
These data are consistent with published findings by Hilgar and 
Hummel, "The androgenic and myogenic evaluation of steroids 
and other compounds -assay 1", in AG Hilgard. DJ, Hummel 

25 (ed.) Endocrine Bioassav Data. Part III )U,S. Dept. HEW NIH 

1964), that 5a-dihydrotestosterone is a highly active androgen, 
whereas 5p-dihydrotestosterone is essentially inactive. 

In the case of the tetrahydro reduced stereoisomers 
of progesterone, 3a-hydroxy-5P-pregan-20-one was a poor fit 
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whereas its stereoisomer 3a-hydroxy-5a-pregnan-20"One was an 
excellent fit. The former molecule is a highly active neurosteroid, 
as reponed by Purdy, et al., /. Med Chem. 33:1572-1581 (1990). 
Comparison of the specific pattern of donor/acceptor linkages of 
3a-hydroxy-5a-pregan-2{)-one with those in the steroid/thyroid 
hormone superfamily demonstrated that the linkage pattern of 3a- 
hydroxy-5a-pregan-20-one is unique, Con5)ounds having such 
unique linkages are predicted to have unique biologic function and 
may be ligands for newly discovered "orphan receptors." 

Example 4 

Correlation of Toxicity and ''Side Effects" with Fit 

into DMA. 

Ligands that fit into more than one site in DNA have 
been observed to have multiple biologic actions. Both desirable 
and undesirable "side effects" should thus be predictable from the 
specific DNA sequence which a given compound fits into as well 
as the manner and relative degree of fit. Examples of molecules 
that have been observed to fit into more than one site in DNA 
include the psychotropics cocaine, morphine, LSD and 
tetrahydrocannabinoids, and certain intercalating antibiotics. For 
example, the monoamine oxidase inhibitor selegiline fits into the 
site in DNA which acconmodates glucose and various oral 
antidiabetic dmgs as reported by Rowland et al., /. Clin, 
Pharmacol 34:80-85(1994). This observation is consistent with 
the finding that selegiline causes hypoglycemia in some patients. 
Another case is the anti-androgen anandron which fits into the site 
in DNA which acconunodates testosterone. Because anandron fits 
into DNA in two orientations, i.e., in a manner sinailar to both 
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androgens and anti-androgens, mixed activity is predicted for diis 
compound. Experimental results indicate that anandron has bodi 
agonist and antagonist activities, as reported by Steinsapir, et al.. 
The Endocrine Society (74th Annual Meeting) 1992:109 (abs. 
228). 

It has also been noted that ligands which cause stress, 
chemical modifications and/or covalent linkages to the DNA when 
fit into a given site fi-equently possess toxicity. Examples include 
certain carcinogens and teratogens, e.g., thalidomide, dioxin, 
arene oxides, aflatoxins and some diethylstilbestrol metabolites. 
Another example is the anti-progestin RU486 which stresses base 
pair hydrogen bonds when inserted fully into the progesterone 
site in DNA. Similar strain is produced by other anti-progestins 
having the same side chains (e.g., llp-phenylamines) raising the 
possibiUty that such features may correlate with abortifacient 
activity attributed to RU486 and related analogs. Observations 
with thahdomide enantiomers indicate that teratogenicity 
associated with this compound may correlate with a stereospecific 
effect on base pairing. An exann5)le of such a toxicity 
pharmacophore, called a toxicophore, is presented in Figures 4E 
and 5E and the corresponding data file is submitted on magnetic 
tape. 

Example 5 

Design and Development of New Drugs, 
Using the principles described above and in the 
examples, new drugs can be designed or existing drugs can be 
redesigned while at the same time Umiting potential undesirable 
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side effects. One example of how an active drag can be designed 
using the technology follows: 

3-Phenylacetylaniino-2,6-piperidinedione (AlO) is a 
modified amino-acid derivative, which was originally isolated 
from fireeze-dried htmMn urine. Despite having low toxicity, 
high concentrations of AlO w»e required to demonstrate 
significant growth inhibitory activity on tumor cells. The focus 
of the following study was to develop more potent analogs. 
Modeling studies demonstrated that AlO was capable of inserting 
into partially unwound double stranded DNA and forming a 
single hydrogen bond between the imino proton of the 
piperidinedione ring and a phosphate oxygen on a single strand. 
It was observed that placing a hydroxyl group at the para position 
of the phenyl group of AlO would enable a second hydrogen bond 
to form thereby substantially enhancing fit as reported by Hendry 
et al., /. Steroid Biochem. Molec. Biol 48:495-505 (1994), the 
teachings of which are hereby incorporated by reference in their 
entirety. The relative fit of AlO (normalized to 100%) and 
various related analogs measured by energy calculations 
demonstrate that the best fitting compound is p-OH-AlO (164%). 
Subsequent synthesis of fliese compounds followed by testing in 
various animal and human tumor cells demonstrated that p-OH- 
AlO was the most active compound and was as much as an order 
of magnitude more active than AlO, as reported by Hendry, et al.. 
Recent Advances in Chemotherapy, Buchner and Rubinstein (eds) 
2498-2499 (1991), Hendry et al., U.S. Patent 5,238,947 which is 
incorporated herein by reference. 
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Example 6 

Comparison of the Drug Design Technology to 
Classical Structure-Activity Methods, 

The drug design technology described here can be 
used in conjunction with quantitative-structure-activity- 
relationship methods (QSAR), e.g., comparative field molecular 
analysis (CoMFA). One value of the approach is that it facilitates 
the orientation of various ligands relative to one another in three 
dimensions. The successful structure-activity relationship found 
for estrogens derived from fit into DNA is described here. If one 
were to attempt to derive such a relationship a priori without first 
knowing the detailed tiiree-dimensional stracture of an 
appropriate macromolecule (e.g., the ligand binding site of a 
receptor or an enzymatic site), chemical intuition would 
necessitate searching for conmaon features that exist in known 
active structures. In the case of the natural hormone estradiol and 
the potent synthetic estrogen trans-diethylstilbestrol, such a 
common feature is the phenoxy group. Alterations of the 
phenoxy group give rise to inactive stractures. If one overlaps 
the resulting three dimensional orientation with that which is 
obtained by optimal docking of these molecules into DNA, a 
different pattern emerges. Thus, one would expect that such 
different orientations when used subsequently to correlate activity 
of otiier molecules would give very different results. In fact, 
analysis of relative fit into DNA in kcal shows diat if the 
orientation of trans-diethylstilbestrol based upon overlapping the 
A ring of estradiol is used to dock the ligand into DNA, a 
relatively poor fit results. In this case, poor activity for trans- 
diethylstilbestrol would be incorrectly predicted. In contrast. 
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using the orientation of DES derived only from the 
stereochemistry of DNA as taught by the present invention, 
increased activity would be correctly predicted. 

Example 7 

Design of an anti-estrogen, para- 
hydroxyphenylacetylaminO'2y6'piperidinedione, a regulator of 
tumor cell growth. 

Within the past few years there has been a growing 
interest in nontoxic, naturally occurring small molecules as 
regulators of tumor cell growth. Examples of recently published 
findings include: regression of mammary carcinomas by a dietary 
monocyclic nonoterpens, limonene; modulation of oncogene 
expression in erythroleukemic cells growth by an endogenous 
product of lipid peroxidation, 4-hydroxynonenal; inhibition of 
malignant cell growth by the endogenous ligand p- 
hydroxyphenylacetate; and induction of tumor cell differentiation 
in premalignant and malignant cells by a circulating component of 
human plasma, phenylacetate. 

Phenylacetate has been shown to reduce levels of the 
myc oncogene which is involved in the development of several 
cancers including breast, brain, prostate, blood, lung and colon. 
Another mechanism by which phenylacetate is thought to be 
effective is by reducing levels of the amino acid glutamine. 
Phenylacetate conjugates witfi circulating glutamine to produce 
die excreted urinary metabolite phenylacetylglutamine (PAG). 
Cancer cells require glutamine for growth and are known to be 
more sensitive to glutamine depletion that normal cells. These 
findings have led to the initiation of Phase I clinical trials with 
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phenylacetate in brain and prostate cancer at the National Cancer 
Institute. 

In the process of screening fractions of freeze dried 
human urine for growth inhibitory activity in human breast 
cancer cells, a dehydration product of PAG was isolated. The 
compound was characterized as 3-phenylacetylamino-2,6- 
piperidinedione by spectroscopic methods and independent 
synthesis and was termed AlO based upon the chromatographic 
fraction from which it was isolated. It has not been conclusively 
determined whether this structure is a circulating compound, 
however, it is similar to phenylacetate in that it lacks toxicity in 
both laboratory animals and humans. The compound has also 
been found to inhibit growth in a variety of cancer cells in vitro, 
as well as human breast cancer transplanted into athymic nude 
mice. Oiemoprevention effects have also been reported. In 
addition to phenylacetate and methyl p-hydroxyphenylacetate, 
there is a number of synthetic compounds with antitumor 
activities that have structural features in common with 3- 
phenylacetylamino-2,6-piperidinedione. Examples of such 
compounds in which analogies to the piperidinedione ring are 
prominent include the aromatase inhibitors aminoglutethimide, 
rogletimide and related analogs which are used in the treatment of 
breast cancer; the alkylating drug PCNU which inhibits tumor 
growth by proposed interaction with DNA; 5-cinnamoyl-6- 
aminouracil derivatives which inhibit tumor growth via putative 
DNA intercalation; amonafide and its congeners which mediate 
topoisomerase n DNA cleavage by intercalation, bis (2,6- 
dioxopiperazine) derivatives, e.g., ICRF-193, which are potent, 
direct inhibitors of mammalian DNA topoisomerase E. 
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Relative high doses of 3-phenylacetylainino-2,6- 
piperidinedione have been generally required to inhibit tumor 
growth both in vitro and in vivo. The goal of this study was to 
identify more active phenylacetylamino-2, 6-piperidinediones 

5 using the technique described above. The technology is based 

upon modeling of the stereospecific fit of molecules into DNA 
and has been recently modified to take advantage of computer 
graphics and energy calculations. Computer modeling was 
followed by the design, synthesis, and in vitro biological testing of 

10 3-phenylacetylamino-2,6-piperidinedione derivatives. The 

molecule predicated to be the most active based upon degree of fit 
in DNA, i.e., the p-hydroxy derivative, was found to be the most 
active antitumor agent in all of the biological assays investigated. 
When tested in MCF-7 (E3) human breast cancer cells, the p- 

15 hydroxy derivative possessed antiestrogenic activity in die range 

of the drug tamoxifen which is currently in clinical use for the 
treatment of breast cancer. 

Example 8 

Three separate molecules designed with modeling 
20 technology are shown in Figure 7. These molecules, termed SGI 

100, SGI 101 and SGI 102 bear similarities and differences to 
each otiier. They all show structural similarities to components of 
both estrogen and progesterone. SGI 100 was designed on the 
basis of its abihty to fit into the site in DNA which accommodates 
25 both estradiol and progesterone. The manner in which it fits 

predicted antagonist activity. When fit to die estrogen 
pharmacophore (Figure 2 and Appendix 1), die acetyl group at 
the 17p position extended out of die pharmacophore and had 
electrostatic repulsion witii dummy atoms (-131 kcal) compared 
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to estrogen (-51 kcal). Accordingly, this high positive energy of 
interaction indicates that SGI 100 acts as an antagonist. The 
binding of SGI 100 to the estrogen receptor is dose dependent and 
approximately 144 times less than estradiol In bioactivity 
experiments, SGI 100 significantly decreased cell growth in MCF- 
7 human breast cancer cells (134,431 cells) at a dose of 10 * M 
when compared to control cells (252,197 cells). The same 
concentration of tamoxifen citrate decreased the number of MCF- 
7 cells to 187,759. Thus, the design of this compound based on the 
pharmacophore approach of this invention predicted a 
demonstrable anti-estrogen bioactivity that was greater than 
tamoxifen. SGI 1(X) binds in a dose dependent manner to the 
progesterone receptor but with 133 to 200 times less affinity. 

SGI 101 was designed on the basis of fitting into 
DNA at the site which acconraiodates estradiol but with opposite 
hydrogen bonding properties which predict estrogen antagonist 
activity, SGI 101 extends beyond the estrogen phannacophore and 
has electrostatic repulsion between the para-nitro group and 
dummy atoms. SGI 101 is the most potent analog designed by the 
technology as measured by growth inhibition of MCF-7 cells. At a 
dose of 10"* M, SGI 101 inhibited cell growth (81,103 cells) 
relative to control (252,197 cells) and was substantially more 
active than the same concentration of tamoxifen (187,759 cells). 

SGI 102 was designed on die basis of its fit into DNA 
at the site which accommodates progesterone. SGI 102 possesses 
an alkyl amino side chain at the 1 ip position which extends out of 
the site between base pairs into the major groove. It has different 
hydrogen bonding properties than progesterone and would extend 
beyond the volume map of the progesterone phannacophore. As 
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such, it is predicted to be an antagonist. SGI 102 was designed 
prospectively, synthesized, and tested in various biological assays. 
SGI 102 binds in a dose dependent manner to the progesterone 
receptor but not as strongly as progesterone or the abortifacient 
antiprogestin RU486. In animal experiments, SGI 102 showed no 
abortifacient activity. However, in experiments using MCF-7 
human breast cancer cells, SGI 102 had equivalent activity to 
RU486 in inhibiting growth. These findings are consistent with 
the predictions made by the modeling technology. 

It should be imderstood diat the foregoing relates 
only to a preferred embodiment of the present invention and that 
numerous modifications or alterations may be made therein 
without departing from the spirit and the scope of the invention as 
set forth in the appended claims. 
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CLAIMS 

I claim: 

1. A method for designing molecules having 
biological activity comprising 

designing compounds which are complementary to a 

double stranded DNA helical structure; 
determining the energy of interaction between the 

compoimd and the DNA in the presence of specific 

biological receptors; and 
selecting compoimds having an equal or lower energy 

of interaction than the natural ligands for the 

specific biological receptors. 

2. The method of Claim 1 wherein the receptor 
binds to a compoimd selected firom the group consisting of 
steroids, thyroid hormones, vitamins, phytohormones, peptide 
hormones, neurotransmitters, prostaglandins, and sugars. 

3. The method of Claim 1 wherein the compound 
has a cyclopentanophenanthrene motif. 

4. The method of Claim 1 wherein the compound 
is complementary to an unwound cavity within a DNA helix. 

5. The method of Claim 1 wherein the compound 
is complementary to a cavity in DNA formed by shifting a base in . 
an unwound site. 
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6. The method of Claim 1 wherein the compound 
is complementary to a cavity formed in DNA by removing a 
nucleotide base. 

5 

7. The method of Claim 1 wherein the compound 
is a receptor antagonist. 

8. The method of Claim 7 wherein the conq)ound 
10 is reproductive hormone antagonist. 

9. A compound having biological activity isolated 
by the process comprising, 

designing compounds which are complementary to a 

double stranded DNA helical structure; 
determining the energy of interaction between the 
compound and the DNA in the presence of specific 
biological receptors; and 
selecting those compounds having an equal or lower 
energy of interaction than the natural ligands for the 
specific biological receptors. 

10. The compoimd of Claim 9 wherein the 
receptor binds to a compound selected from the group consisting 

25 of steroids, thyroid hormones, vitamins, phytohormones, peptide 

hormones, neurotransmitters, prostaglandins, and sugars. 

11. The compound of Claim 9 wherein the 
compound has a cyclopentanophenandarene motif. 



15 



20 
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12. The compound of Claim 9 wherein the 
compound is complementary to an unwound cavity within a DNA 
helix. 

13. The compound of Claim 9 wherein the 
compound is complementary to a cavity in DNA formed by 
shifting a base in an unwound site. 

14. The compound of Claim 9 wherein the 
confound is complementary to a cavity formed in DNA by 
removing a base. 

15. The compound of Claim 9 having estrogenic 

activity. 

16. The compound of Qaim 9 wherein the 
compound is a receptor antagonist. 

17. The compound of Qaim 16 wherein the 
compound is a reproductive hormone antagonist. 

18. A pharmacophore comprising, 

a three-dimensional array of points defining a specific 
shape and volume, wherein said three-dimensional 
array of points is the aggregate average shape of a 
plurality of molecules when optimally fit into an 
unwound DNA site, said molecules possessing the 
same or similar biological activity. 
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19. A biologically active compound designed using 
a pharmacophore of Claim 18, wherein tiie template identifies 
areas of structure and charge responsible for fit into DNA, 

20. A method for creating a pharmacophore 

comprising the steps of : 

determining tiie optimal fit of a plurality of 
compoimds having the same or similar biological 
activity into nucleic acid sequences such that the 
lowest energy of interaction and best steric fit are 
obtained; 

aligning the compounds relative to the heteroatoms on 

the nucleic acids; and 
defining a diree dimensional shape representing die 

aggregate average shape resulting from said 

alignment of the compounds. 

21. The metiiod of claim 20, wherein the nucleic 
acid sequence is selected from die group consisting of 
deoxyribonucleic acid, double stranded deoxyribonucleic add, 
ribonucleic acid, complexes of deoxyribonucleic acid and 
ribonucleic acid, and apurinic and apyrimidinic sites. 

22. The metiiod of claim 20, wherein die nucleic 
acid is deoxyribonucleic acid. 
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23. A method of screening a molecule for 
bioactivity or toxicity wherein the molecule is examined for the 
ability to fit into a pharmacophore comprising the steps of: 

measuring the energy of interaction when the molecule 
is fitted into the pharmacophore; and 

comparing the energy of interaction of the molecule to 
a predetermined energy of activation that correlates 
to a biological activity. 

24. The method of claim 23 wherein the 
pharmacophore is selected from the group consisting of an 
estrogen pharmacophore, an antiestrogen pharmacophore, an 
androgen pharmacophore, a thyroid hormone pharmacophore, 
and toxicophore. 

25. A method for screening a molecule for toxicity 
wherein the molecule is examined for the ability to fit into a 
toxicity pharmacophore comprising the steps of 

measuring the energy of interation when the molecule 
is fitted into the pharmacophctt^; and 

comparing the energy of interaction of the molecule to 
a predetermined energy of acdvation that correlates 
to a toxic activity. 
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26. A method of designing a molecule with a 
desired biological activity comprising the step of determining the 
structure of a molecule which will have a minimal predetermined 
energy of activation when fitted into a pharmacophore with a 
5 desired activity. 
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Relative Fit To Tlie Estrogen Pharmacophore 



FIGURE 6 
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